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NOMENCLATURE 


a  time  at  which  the  applied  pressure  impulse  ends 

times  that  define  the  constant  pressure  r  quare  wave  segments 

A  [VgjW-YbJXfc 

hj  weighting  coefficients  determined  from  /  (/) 

B.  7 b, 

D  shock  wave  propagation  speed 

dlc  shock  penetration  depth  in  a  deep  snow  cover 

/  (t)  function  describing  the  shape  of  a  variable  pressure  impulse 

H,  dH  momentum  per  unit  area  and  its  derivative 
Hp  total  momentum  per  unit  area  applied  to  the  snow 

Hp  momentum  per  unit  area  applied  to  the  snow  in  the  ith  square  wave 
1  segment 

Ht  momentu  m  per  unit  area  in  the  snow 

IQ  instantaneously  applied  pressure  impulse 

w,  dm  mass  per  unit  area  and  its  derivative 

P  shock  wave  pressure 

PQ  maximum  pressure  amplitude  of  the  shock 

P(0  shock  pressure  amplitude  as  a  function  of  time 

r  current  position  of  the  inner  radius  of  the  cavity  surface 

R  Bhock  wave  propagation  radius 

Rq  initial  radius  of  a  cavity  surface  in  the  snow 

t,  dt  time  and  its  derivative 

tf  time  required  for  the  shock  to  reach  Xf 

tQ  time  at  which  the  pressure  impulse  is  applied 

U{)  distance  snow  surface  has  moved 

V  snow  particle  velocity 

V{  snow  particle  velocity  in  the  ith  square  wave  segment 

shock  front  position  at  t  =  a 
X^  shock  front  position  at  t  = 

Xf  shock  front  position;  Xf  =  Xf-X;)  for  XQ  =  0 

X,  current  snow  surface  position 

XQ  initial  snow  surface  position 

a  attenuation  coefficient 

P  relative  snow  compaction  (1  -  P0/ Pf) 

Y  P0P/Pf(l-P> 

e  R-Ro 

p(  compacted  snow  density 

p„  Initial  snow  density 
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Estimates  of  Shock  Wave  Attenuation  in  Snow 


JEROME  B.  JOHNSON 


INTRODUCTION 

Shock  wave  attenuation  in  snow  is  affected  by  the  amplitude,  geometry  and 
duration  of  a  shock,  as  well  as  the  mechanical  properties  of  the  snow.  Accordingly, 
these  effects  must  be  taken  into  account  in  any  description  of  shock  attenuation.  The 
mechanical  properties  of  snow  for  high-strain-rate,  large-amplitude  shocks  are  not 
now  well  understood.  Consequently,  attempts  to  estimate  shock  attenuation  in 
snow  have  relied  on  field  measurements  of  pressure  attenuation  from  explosions  in 
snow  and  on  theoretical  constitutive  descriptions. 

Mellor  (1977)  used  the  results  from  field  measurements  of  explosive  detonations 
in  snow  to  estimate  the  attenuating  properties  of  snow  for  spherical  shock  waves.  He 
estimated  that  the  attenuation  of  shock  pressure,  combining  geometric  spreading 
and  internal  dissipation,  near  an  explosion  decays  as  /T4  (R  being  the  propagation 
radius)  and  close  to  the  elastic  limit  it  is  approximately  an  R "3  decay.  Brown  (1980, 
1981,  1983)  appears  to  have  made  the  first  and  only  attempts  to  estimate  shock 
attenuation  for  plane  waves  in  snow.  He  used  two  theoretical  volumetric  constitu¬ 
tive  laws,  one  for  medium  to  high  density  snow  and  one  for  low  density  snow,  to 
estimate  shock  attenuation.  His  results  indicate  that  plane  shock  waves  can  attenu¬ 
ate  by  more  than  80  to  90%  after  propagating  through  only  0.04  to  0.06  m  of  snow. 
These  are  large  attenuations  that  are,  at  their  maximum,  about  an  (Xf-Xu)-1,2  decay 
([Xj-XQ]  is  the  propagation  distance). 

Brown  did  not  discuss  the  importance  of  an  applied  shock's  amplitude,  geom¬ 
etry  and  duration  on  attenuation  features.  This  paper  will  examine  the  importance 
of  the  pressure-time  profile  of  an  applied  pressure  impulse  on  shock  attenuation  in 
snow.  Snow  will  be  represented  by  a  simple  mechanical  model  i  n  which  attenuation 
occurs  through  momentum  transfer  from  the  applied  shock  !  >  the  snow. 


MOMENTUM  MODEL 

The  momentum  model,  also  known  as  the  "snowplow"  model,  was  used  in  initial 
attempts  at  developing  constitutive  relations  to  analyze  the  dynamic  behavior  of 
porous  materials  (Herrmann  1971),  The  porous  material  is  assumed  to  compact  to 
its  final  density  at  a  negligible  stress.  After  compaction,  the  material  is  assumed  to 
be  incompressible  (an  ideal  locking  material).  The  change  in  momentum  caused  by 
a  stress  wave  (the  pressure  impulse)  is  assumed  to  be  spread  uniformly  over  the 
volume  of  the  material  behind  an  advancing  shock  wave,  This  means  that  the  stress 
wave  is  lengthened  in  time  and  reduced  in  amplitude  as  more  of  the  material  is 
compacted  by  the  propagating  shock.  Hence,  attenuation  is  caused  by  momentum 
spreading,  and  losses  attributable  to  plastic  deformation,  fracturing  and  release 
waves  are  not  considered. 


INSTANTANEOUSLY  APPLIED  IMPULSE 


Consider  a  pressure  impulse  tf0)  applied  normal  to  the  plane  surface  of  snow,  and 
assume  that  snow  is  an  ideal  locking  material.  Snow  next  to  the  surface  will  be 
immedia  tely  compacted  to  its  final  density.  Since  the  compacted  snow  is  rigid,  it  will 
move  at  a  uniform  pressure  and  particle  velocity  after  the  pressure  impulse  is  ap¬ 
plied.  The  stress  wave  will  propagate  into  the  snow  as  a  compaction  shock,  moving 
with  a  velocity  D  at  a  pressure  P  and  particle  velocity  V.  At  the  shock  front  these 
parameters  are  related  by  the  Rank!  ne-Hugoniot  jump  conditions  for  the  conserva¬ 
tion  of  mass  and  momentum  across  the  shock  front  (Kolsky  1963) 

PoDapf  [D-V]  (1) 

and 

P  »  P0DV  (2) 

where  p0  is  the  initial  density  of  the  snow  and  pf  is  the  compacted  snow  density.  The 
compacted  snow  density  pf  is  usually  determined  by  experiment,  but  is  estimated 
for  this  paper  because  of  the  lack  of  suitable  da  ta.  In  using  eq  1  and  2,  the  contribution 
of  shock  heating  in  the  snow  is  neglected. 


p* — u0— ^ 


x0  X|  X, 

Figure  1 .  Deformation  geometry  for  plane  shock  mve  propa¬ 
gation  in  snow. 


Figure  1  shows  the  deformation  geometry  in  one  dimension  for  a  pressure  im¬ 
pulse  applied  at  the  initial  snow  surface  XQ.  The  location  of  the  surface  of  the  snow 
and  the  shock  front  at  some  time  after  the  application  of  the  pressure  impulse  are 
and  Xf  respectively.  During  the  time  that  the  shock  front  has  traveled  to  X(,  the  snow 
surf  ace  has  moved  a  distance  UQ.  UQ  is  a  function  of  time  or,  alternatively,  a  function 
of  shock  propagation  position  Xf .  At  any  time  the  location  of  the  snow  surface  is 

X|  =X0  +  U0.  (3) 

The  displacement  of  the  snow  surface  can  be  calculated  by  integrating  the  parti¬ 
cle  velocity  over  time 

tf 

U0  =  Vdt  (4) 

•  to 

where  tQ  is  the  time  that  the  shock  is  applied  and  ff  denotes  the  time  required  for  the 
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shock  to  reach  X(..  Knowing  the  shock  pressure  attenuation  with  distance  is  of  more 

practical  interest  than  following  the  pressure  with  time.  Thus,  the  snow  surface 
displacement  will  be  reformulated  in  terms  of  shock  position  rather  than  time.  By  use 
of  the  fact  that 


D  =  dlL 
dt 

eq  4  may  be  rewritten  as 


Uo 


From  eq  1 


(5) 


(6) 


(7) 


where  p  describes  the  relative  snow  compaction.  Substituting  eq  7  into  eq  6  and 
solving  gives  the  position  of  the  snow  surface  at  time  tf 

X|  =(l-p)  X0  +pXf.  (8) 

The  mass  and  momentum  per  unit  area  of  the  snow  per  unit  length  contained  in  the 
length  (Xf  -  Xj)  between  the  current  snow  surface  position  X(  and  shock  front 
position  Xf  are 

dm  =  Pf  dx  (9) 


dH  -  P,V  dx  .  (10) 

Since  the  compacted  snow  is  assumed  to  be  rigid,  V  is  constant  through  the  region 
X(  to  Xf  and  is  equal  to  the  particle  velocity  at  X{,  i.e.  V  =  /(X,).  Consequently,  the 
momentum  per  unit  area  in  the  snow  at  the  time  the  shock  has  reached  p(  is 


H. 


fXf 

PfVdx  =>Vpf  {Xf  -X,). 

Jx, 


(11) 


Substituting  in  for  X)  from  eq  8  gives 


H,  =PfV(l-p)(Xf  -X„)  . 


(12) 


The  momentum  per  unit  area  in  the  snow  Hb  must  be  equal  to  the  momentum  per 
unit  area  applied  to  the  snow  from  the  instantaneously  applied  pressure  impulse 
that  is  given  by 


Hp  -  lo- 

Equating  and  and  solving  for  V  gives 


(13) 
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Equations  7  and  14  can  be  used  in  eq  2  to  determine  the  pressure  at  the  shock  front 


P(Xf)  = 


P0P('Xf  “Xq f 


(15) 


Figure  2  snows  the  shock  front  pressure  for  snow  subjected  to  an  instantaneous 
pressure  impulse  of  IQ  =  100  Pa  s.  The  pressure  attenuates  as  (Xf-XQ)2,  with  pro¬ 
portionality  constants  that  depend  on  impulse  magnitude,  initial  snow  density  and 
the  relative  snow  compaction  p.  Figure  2a  shows  that  snow  with  a  higher  initial 
density  supports  less  pressure  at  the  same  propagation  distance  as  compared  to 
lower  density  snow, 

Equation  15  can  be  rewritten  to  show  shock  pressure  as  a  function  of  the  total 
mass  of  compacted  snow 

P  (Xf)  m  llh  (15a) 

pm2 

where 


m  •  p0(Xf~XQ). 


(15b) 


a,  Initial  demit  ies  of  100, 250  and  400  kg/m3  and  final 
density  of  900  kg/m3. 


b.  Initial  density  of  200  kg/m3  and  final  densities  of 
400, 600  and  900  kg/m3. 


Figure  2.  Pressure  decay  with  distance  for  an  instantaneously  applied  plane  shock  wave  pressure  impulse  (100  Pa 
s),  predicted  by  the  snowplow  model. 
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Figure  3.  Pressure  decay  as  a 
function  of  snow  mass  com¬ 
pacted  by  the  propagating 
plane  shock  wave ,  predicted 
by  the  snowplow  model  for  the 
conditions  given  in  Figure  2a. 


The  results  given  in  Figure  2a  are  recalculated  using  eq  15a  and  shown  in  Figure  3. 
Figure  3  shows  that  snow  with  an  initially  lower  density  will  support  a  lower 
pressure  magnitude  than  will  a  snow  with  a  higher  initial  density,  given  the  same 
total  compacted  mass.  This  finding  leads  to  the  common  observation  that  lower 
density  materials  are  better  attenuators  than  higher  density  materials  because  less 
mass  is  needed  to  produce  a  given  attenuation  than  for  a  porous  material  of  higher 
density. 


CONSTANT  PRESSURE  IMPULSE 
OF  FINITE  DURATION 

The  solution  given  by  eq  15  is  not  a  very  satisfying  way  of  evaluating  shock 
attenuation  since  it  assumes  instantaneous  application  of  the  pressure  impulse.  This 
is  unrealistic,  as  it  implies  infinite  pressure  at  x  =  X0,  and  does  not  show  how 
pressure  pulses  with  the  same  total  impulse  but  different  amplitudes  and  durations 
are  attenuated.  An  analytical  solution  can  be  derived  for  a  constant  pressure  impulse 
of  finite  time  duration,  given  by 

Pit)  »  P0  0  £  t  £  a 

(16) 

P(f)  »  0  a  <  t 

where  a  is  the  time  duration  of  the  applied  pressure  pulse.  The  momentum  per  unit 
area  caused  by  the  pressure  impulse  is 

Hp  =  f  P0  dt  for  0  £  f  £  <?  (17) 

Jo 


and 
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(a 


H  p- 


P0  iff  tor  a  <  t 


(18) 


Equations  17  and  18  show  that  the  momentum  during  application  of  the  pressure 
impulse  will  vary  with  time,  while  the  momentum  after  the  impulse  has  been 
applied  is  constant.  Therefore,  separate  solutions  for  shock  pressure  attenuation  are 
needed  during  the  time  period  of  pulse  application  and  for  the  time  period  after  the 
pulse  has  been  applied.  Equations  17  and  18  can  be  transformed  into  spatial 
coordinates  by  use  of  the  results  of  eq  5  and  7,  giving 


Hp 


for  0  <  Xf  <  <&-  . 

P 


(19) 


The  limits  of  integration  are  0  to  *  a  V  /  jJ.  Equa  ting  momentum  in  the  snow  (eq  1 2) 
to  the  pressure  impulse  momentum  (eq  19)  and  differentiating  gives 


-fofi-r  »  V  -f-  [V  X f] .  (20) 

Pf(l-P)  d* 

Also 


<21) 

Substituting  the  identity  in  eq  21  into  20  and  integrating  gives  the  particle  velocity 


for  0  SXf  <;  Xa. 


(22) 


Hence,  during  the  period  of  pressure  impulse  application,  the  particle  velocity  is 
constant  and  there  is  no  attenuation  of  the  shock  pressure.  The  distance  that  the 
shock  has  traveled  at  t  =  a  is  given  b> 


Xa  a  Da 


Ppti^ _ 

.PPf(l-P). 


1/2 


(23) 


After  the  pressure  impulse  momentum  has  been  applied  to  the  snow,  the 
relationship  between  the  pressure  impulse  momentum  and  snow  momentum  is 


Port  =PfV(l-p)Xf  for  Xa  <  X( . 


(24) 


Solving  for  the  particle  velocity  and  using  the  definition  of  Xfl  in  eq  23  gives 


V  = 


PoP  l1/2Xa 

pf(i-P)J  V 


(25) 


The  shock  pressure  as  a  function  of  distance  can  now  be  determined  from  eq  2, 7, 22 
and  25,  giving 
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P(X|)=P0  for  0  £  Xf  £  X, 


(26a) 


P(X()  -  forXa  <  Xf.  (26b) 

Figure  4  (square  wave)  shows  the  pressure  as  a  function  of  distance  for  an  initial 
snow  density  of 350  kg/  m3,  |3 «  0.3  and  a  pressure  impulse  of 605  Pa  s,  The  shock  does 
not  begin  attenuating  until  after  all  of  the  pressure  impulse  momentum  has  been 
applied  to  the  snow  surface,  and  then  it  attenuates  as  X"?  Comparison  of  the  square 
wave  pressure  impulse  to  an  instantaneous  pressure  impulse  (impulse)  shows  that 
the  pressure  is  low  er  for  the  square  wave  during  its  duration  of  application  than  for 
an  instantaneously  applied  pressure  impulse.  After  the  square  wave  momentum 
has  been  completely  transferred  to  the  snow,  both  the  square  wave  and  instantane¬ 
ously  applied  impulse  pressures  are  the  same.  Since  explosives  arid  other  sources  of 
shocks  are  applied  over  a  finite  time,  this  result  shows  the  importance  of  modeling 
the  pressure-time  profile  of  an  applied  pressure  impulse  for  predicting  attenuation. 


Distance  (m) 

Figure  4.  Comparison  of  pressure  decay  with  distance  for  plane 
shock  waves,  each  having  a  total  applied  pressure  impulse  of 
605  Pas,  predicted  by  the  snowplow  model,  Thepressure  impulse 
was  applied  instantaneously  for  impulse,  with  a  square  wave  of  0,1 7 
ms  duration  for  square  wave,  and  an  exponential  pulse  of  0.1 6  ms 
duration  for  exponential, 


VARIABLE  PRESSURE  IMPULSE 
OF  FINITE  DURATION 

Simulating  a  realistic  applied  pressure  impulse  requires  a  function  that  can 
represent  variable  pressure  impulses  of  finite  duration.  These  variable  pressure 
impulses  result  in  nonlinear  differential  equations  that  do  not  have  closed  form 
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solutions  and  whose  numerical  solutions  can  be  unstable.  In  addition,  the  form  of 
the  differential  equations  depend  on  the  particular  variable  pressure  impulse  ap¬ 
plied  to  the  snow  surface.  To  reduce  these  computational  difficulties,  the  variable 
pressure  impulses  were  approximated  by  a  sequence  of  square  waves.  The  accuracy 
of  such  an  approximation  depends  on  the  time  duration  specified  for  the  square 
wave  segments  that  constitute  the  applied  pressure  impulse.  Thus,  instead  of  ex¬ 
actly  formulating  the  problem  and  obtaining  a  numerical  solution,  the  pressure 
impulse  is  approximated  and  solved  analytically  within  the  sequential  time  steps. 
The  variable  pressure  impulse  of  finite  duration  beginning  at  t  -  0  is  approximated 
by 


P(t)  *P0bi 


bo  =  1 

t,»0 


for  0  £  t  <  a  i 
for  a  i  £  t  <  a  2 

for  a  n_i  £  t  <  a  n 
for  a  n  £  t . 


(27) 


The  shape  of  the  vanable  pressure  impulse  is  fit),  where  the  b,  weighting  coefficients 
are  determined  at  the  midpoint  between  the  beginning  time  and  ending  time  for 
each  square  wave  segment  fl(aM  +  fl])/2].  The  aj  are  the  times  that  define  the  square 
wave  segments  where  i  =l,2,3,...,n-l. 

Solving  for  the  pressure  in  snow  as  a  function  of  distance,  using  eq  27,  requires 
that  the  problem  be  formulated  for  three  conditions.  The  first  is  for  the  time,  or 
spatial  increment,  of  the  first  square  wave  segment.  The  next  solu  tion  must  account 
for  the  momentum  that  has  been  applied  to  the  snow  at  the  time,  or  spatial  position, 
of  interest.  Finally,  a  solution  must  be  found  after  all  of  the  pressure  impulse 
momentum  has  been  applied  to  the  snow.  These  conditions  can  be  expressed  as 


(X( 

Pobo 

.0 


JL  dx 

v0(x) 


PfV0(Xf)(l-p)Xf  for  0  ^  Xf  <  Xa) 


(28a) 


f* 

Mp„  + 

JX a, 


-  PfV,(Xf)(l-P)  Xf 


for  X8)  £  Xf  <  Xflj 


n-2 

I  Hp  + 
l«0  1 


(Xf 


Vn-ir^r^  «  PfVVl  (Xf)(i -P)Xf 


(28b) 


for  Xfln  l  S  Xf  <  Xfln 


n-l 

I  Hn  =  Pf  V„(Xf)(l  -P)Xf  for  Xa  S  Xf .  (28c) 

i-0  1 
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There  are  /;+l  equations  to  solve  before  the  complete  particle  velocity  solution  can 

be  found.  Each  equation  depends  on  the  solutions  of  the  previous  segments. 

Equation  28a  is  a  constant  pressure  impulse,  with  the  solution  given  by  eq  22, 23 
and  26a 


Vo(Xf)  ■  y  forOSXf<Xai  (30) 


where  y  =  ,-Hfl.lL  , 
pfM 

Equations  28b  can  be  solved  sequentially  to  give  the  solution  for  the  spatial  position 
of  interest.  These  solutions  can  be  given  as  a  set  of  recursion  relations  for  position, 
particle  velocity,  and  time 


Vf(Xf)  (32) 

xf 

where 

i  =  1,2, . . . ,  n-1  for  Xal  £  Xf  <  Xa,  +1 
and 

*1  “[^(Xo^-YfcJXa2, 

B|  as  yb\  . 

n-i 

XHPi 

Vn(Xf)  =  — — - —  for  Xa_  <  Xf  (33) 

V  U  Pf(l-p)Xf  n  1 

where 

n-1 

,5,  Hp>  ■  hp 

is  the  total  momentum  per  unit  area  applied  to  the  snow.  Equations  29  and  30  are 
used  to  start  the  solution  and  the  recursion  relations  are  used  to  follow  the  pro¬ 
gression  of  the  shock. 

Figure  4  (exponential)  shows  the  results  of  using  eq  29  through  33, 2  and  7  to 
calculate  the  pressure  as  a  function  of  shock  wave  propagation  distance  for  an  ex¬ 
ponentially  decaying  pressure  impulse  (total  pressure  impulse  =  605  Pa  s) 

P  =  P0e-«'  (34) 

where  a  is  a  decay  constant  and  P0  is  the  initial  pressure  and  f(t)  =  e-a'.  C  jmparison 
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of  the  three  different  applied  pressure  impulses  (each  with  the  same  total  momen¬ 
tum)  in  Figure  4  shows  that  shock  attenuation  can  be  markedly  different  during  the 
period  of  pressure  impulse  application.  The  instantaneously  applied  pressure 
impulse  applies  all  of  its  momentum  at  once  resulting  in  infinite  initial  pressure  that 
immediately  begins  attenuating  as  Xj2,  where  XQ  =  0.  The  square  wave  maintains 
a  constant  pressure  equal  to  the  applied  pressure  until  all  of  its  momentum  has  been 
transferred  to  the  snow  and  then  begins  to  attenuate.  Finally,  the  exponential  pulse 
gradually  attenuates  while  its  momentum  is  being  transferred  to  the  snow,  asymp¬ 
totically  approaching  an  Xf2  decay.  Once  a  pressure  impulse  has  been  applied  to 
the  snow,  pressure  attenuation  is  controlled  by  the  mechanical  properties  of  the 
snow  and  the  magnitude  of  total  pressure  impulse. 


ATTENUATION  OF  CYLINDRICAL  AND 
SPHERICAL  GEOMETRY  SHOCK  WAVES 

Geometric  spreading  can  greatly  Increase  shock  attenuation  in  snow.  Torvik 
(1971)  developed  equations  describing  pressure  decay  for  cylindrical  and  spherical 
geometry  shock  waves  using  the  snowplow  model  and  assuming  instantaneous 
application  of  a  pressure  impulse.  For  a  cylindrically  spreading  shock,  the  pressure 
is  given  by 

P(R)  p  (35) 

PoM2(R-r)2 

/■  =  [pK2  +  (l-p)Hg],/Z.  (36) 


The  pressure  for  a  spherically  spreading  shock  is 


P(R)  m  ; 

PoP  R*(R-r )2 


(37) 


Figure  5,  Deformation  geometry  for  diverging 
shock  waves. 


r  =[pR:'  +  (l-p)Rg]  .  (38) 

/0  is  the  instantaneously  applied  pressure  impulse,  RQ  is  the  Mittal  radius  of  the  cavity 
surface  in  the  snow  on  which  the  pressure  impulse  is  applied,  and  r  is  the  location 
of  the  inner  radius  of  the  cavity  at  some  time  after  the  application  of  the  pressure 
impulse  (Fig.  5). 

When  R  «*  RQ  +  e,  where  e  «  RQ,  a  cylindrical  geometry  shock  decays  as  e"2,  in¬ 
creasing  to  a  decay  of  R-4  for  e  »  R0  (Fig.  6a).  A  spherical  geometry  shock  decays 


Ro-l.Om 


R0-  0.01  \pflM  3S0  kg/m 


R„.  0.01 


R0»  1.0  m 


P0  -3S0k( 

(it  at 


a.  Cylindrical shockwaves  xollh  initial  radii  of 0.01  and  1.0  b,  Spherical  shock  waves  ultli  Initial  radii  of  0.01  and  1,0 

in,  predicted  by  the  snowplow  model,  nt  predicted  by  the  snowplow  model, 


p0-  360  kfl/m 
.  Pi  ■  6(  0 
-  0.3 

200  xSk 

600  ^5Sk  R„-  1 


p0«  350  kg/m 
pt  -  500 

3-0.3 

R„-  1  i 


c,  Cylindrical  shock  waves  in  snow  with  initial  densities  of  d.  Spherical  shock  waves  in  snow  with  initial  densities  of 

200  and  350  kg/tn3  and  initial  radii  of  1  m.  200  and  350  kg/m 3  and  initial  radii  of  1  m, 

Figure  6.  Comparison  of  pressure  decay  with  distance. 


as  e-2  for  e  «  R{y  increasing  to  R"6  for  e  »  R()  (Fig.  6b).  These  findings  show  that 
pressure  decay,  including  the  effects  of  geometric  spreading,  vary  significantly  de¬ 
pending  on  the  initial  radius  over  which  the  pressure  impulse  is  applied  and  the 
distance  from  the  initial  radius.  Figures  6c  and  d  show  the  effect  of  changes  in  initial 
snow  density  and  (i  on  pressure  attenuation  for  cylindrical  and  spherical  waves. 


PRESSURE  DECAY  OF  A 
LINE  CHARGE  ON  SNOW 


Line  charges  are  used  to  clear  minefields  and  to  remove  snow  cornices  from 
mountain  ridge  tops.  Therefore,  it  is  of  practical  interest  to  know  the  extent  that  snow 
reduces  the  effectiveness  of  line  charge  detonations.  For  a  line  charge  resting  on  or 
above  a  snow  cover,  the  cylindrical  pressure  wave  that  hits  the  Bnow  surface  has 
propagated  primarily  through  the  air. 

Ford  (1986)  has  measured  air-blast  pressure,  positive  phase  duration  and  total 
pressure  impulse  at  the  ground  as  a  function  of  lateral  distance  from  the  axis  of  a  line 
charge  explosion  (Table  1 ),  The  positive  phase  duration  is  defined  as  the  duration  of 
the  compressive  (positive)  shock  pressure.  Although  air-blast  pressures  over  snow 
would  differ,  it  is  assumed,  for  this  example,  that  the  air-blast  pressures  given  in 
Table  1  would  also  occur  over  snow,  Furthermore,  the  snow  is  assumed  to  be  deep 
and  shock  pressures  are  determined  only  for  small  propagation  depths  into  the 
snow  so  that  e  « 1?0 ,  where  R0  is  the  distance  from  the  line  charge  to  the  point  of 
interest  on  the  snow  surface,  W ith  these  conditions,  the  snow-ground  boundary  can 
be  neglected  and  shock  pressure  in  the  snow  can  be  estimated  from  the  equations  for 
a  one-dimensional  variable  pressure  impulse  of  finite  duration  (eq  29-33, 2  and  7) 
and  an  exponentially  decaying  pressure  pulse  (eq  34),  Figure  7  shows  the  shock 
pressure  as  a  function  of  distance  from  the  line  charge  at  the  snow  surface  (using 

the  pressure  impulse  data  from  Table  1).  Figure  7  also  shows  the  calculated  pres¬ 
sures,  assuming  p0  ®  350  kg/m3  and  p,  ■  500  kg/m3 ,  as  a  function  of  d,c  for  shock 
propagation  depths  into  the  snowof  0.06, 0,1, 0,2  and  0.3  m.  It  is  evident  that  air-blast 
shock  pressures,  for  a  given  penetration  depth  into  the  snow,  do  not  attenuate  the 
same  amount  for  a  given  shock  penetration  depth  at  the  different  rfLC.  Shock 
attenuation  (that  is,  pressure  reduction  after  a  given  shock  propagation  depth)  is 
greatest 
23.97  m) 

(positive  phase  duration)  increases  with 
d,Q  (Table  1).  The  increasing  positive 

phase  duration  allows  the  shock  to  propa*  Table  1.  Air-blast  measurements— 

gate  through  a  greater  depth  of  snow  shot  2  (after  Ford  1986). 
before  starling  the  e“2  decay.  These  find- 


near  the  line  charge  (rfLC  a  0.96  m)  and  least  at  the  farthest  distance  Wlcb 
.  This  occurs  because  the  duration  of  the  applied  air-blast  pressure  impulse 


ings  suggest  that,  for  shallow  shock  pene¬ 
tration  depth  into  a  snow  cover,  line 
charge  pressure  reduction  depends  on 
the  positive  phase  duration  in  addition 
to  the  total  pressure  impulse,  Also,  shock 
pressure  attenuation  In  a  snow  cover  can 
be  overcome,  somewhat,  by  increasing 
the  positive  phase  duration  of  a  given 
pressure  impulse, 
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Figure  7,  Pressure  decay  with  distance  from  the  line  charge 
for  different  shock  propagation  depths  into  the  snow ,  pre¬ 
dicted  by  the  snowplow  model,  The  air-blast  curve  is  drawn 
from  data  given  in  Table  1 . 


DISCUSSION  AND  CONCLUSIONS 

The  snowplow  model  has  a  limited  ability  to  describe  shock  propagation  and  at¬ 
tenuation  in  snow  because  of  the  simplifying  underlying  assumptions.  The  model 
is,  however,  capable  of  giving  conservative  estimates  of  shock  attenuation  in  snow 
and  of  illustrating  some  of  the  Important  features,  Actual  shock  pressure  attenu¬ 
ation  will  be  less  than  that  predicted  by  the  snowplow  model  using  only  one 
compaction  step  to  the  final  snow  density.  Better  estimates  of  shock  pressure 
attenuation  in  snow  may  be  possible  using  the  snowplow  model  and  a  more  rea¬ 
sonable  compaction  path  for  snow. 

The  snowplow  model  predicts  that  snow  can  cause  shock  pressures  to  decay  as 
a  function  of  Xj2  for  plane  waves,  fT4  for  cylindrical  waves  and  R~*  for  spherical 
waves.  Brown's  (1980)  pressure  decay  estimate  of  about  Xj"1^  for  plane  waveB  is 
different  from  thedecay  predicted  by  the  snowplow  model,  Agreementbetween  the 
snowplow  calculations  and  Brown's  estimates  might  improve  if  n  realistic  snow 
compaction  model  were  used  in  the  snowplow  model. 

Mellor  (1977)  estimates  that  spherical  geometry  shocks,  where  r  »  R(),  attenuate 
as  X  '3  to  P."4  as  compared  to  the  snowplow  model  prediction  of  R~tl,  It  may  be  that 
the  snowplow  model  predictions  are  grossly  off  because  of  the  simple  snow 
compaction  path  used  in  the  model  or  that  the  field  measurements  used  by  Mellor 
to  make  his  estimates  are  not  the  rerult  of  shod'  propagation  in  snow.  The  pressure 
sensors  used  in  many  of  the  field  tests  were  not  able  to  survive  near  source  pressures 
and  were,  consequently,  located  well  outside  the  zone  of  extreme  shock  induced 
snow  compaction.  Outside  the  zone  of  snow  compaction,  pressure  decay  is  due  to 
viscous  dissipation  and  geometric  spreading,  which  produce  much  less  pressure 
attenuation  than  the  pore  collapse  mechanisms  described  by  the  snowplow  model. 
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This  implies  that  Mellor's  estimates  of  pressure  decay  in  snow  for  spherical  shocks 
may  be  too  low,  but  a  final  conclusion  is  not  possible  until  better  snow  compaction 

paths  are  used  in  the  model  calculations. 

Estimating  shock  attenuation  requires  some  thoughtful  application  of  the  results 
of  this  study,  Explosive-induced  shocks  in  snow  often  result  from  an  explosive 
charge  detonation  on  the  snow  or  in  the  air  above  the  snow.  In  this  situation  the 
shock  will  be  transmitted  into  the  snow  by  an  air-blast  wave  propagating  over  the 
snow  surface,  If  the  initial  radius  for  the  air-blast  pressure  wave  striking  the  snow 
is  assumed  to  be  the  distance  from  the  charge  source  to  the  position  at  which  the 
pressure  is  transmitted  into  the  snow,  then  in  most  caseB  RQ  will  be  much  greater 
than  the  snow  depth  (i.e„  e  «  R0),  Attenuation  from  the  transmission  position 
through  the  snow  will  proceed  as  e"2  rather  than  as  ft-5  for  spherical  waves  or  R~* 
for  cylindrical  waves. 

Pressure  decay  of  shock  waves  in  snow  can  be  delayed  by  increasing  the  positive 
phase  duration  of  the  applied  pressure  impulse. 
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